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Abstract. There is a major need to better constrain nuclear parton densities in 
order to provide reliable perturbative QCD predictions at the LHC as well as to probe 
possible non-linear evolution at small values of x. In these proceedings, we discuss 
how the production of prompt photons at large rapidity in p-p and cZ-Au collisions at 
RHIC (y/s NN = 200 GeV) is sensitive to the nuclear modifications of gluon distributions 
at x = O (10 -3 ) and at rather low scales, Q 2 ~ 10 GeV 2 . The nuclear production 
ratio, R dA — a(d + A — * 7 + X )/ (2 A x a(p + p — ► 7 + X )), is computed for isolated 
prompt photons at NLO using the nDSg nuclear parton densities, in order to assess 
the visibility of the signal. We also emphasise that the expected counting rates in a 
year of running at RHIC are large, indicating that R dA could be measured with a high 
statistical accuracy. 



1. Gluons from photons 

An important experimental and theoretical effort has been pursued over the last fifteen 
years aiming at constraining the gluon density in a proton, G P . The current precision 
on G P (x,Q 2 ) has reached a few percent level at not too large nor not too small values 
of Bjorken-x , 10~ 4 < x < 10 _1 , and at large virtualities, Q 2 > 10 GeV 2 . Lacking a 
similar facility such as HERA using heavy ion beams, the gluon density in a nucleus, 
G A (x, Q 2 ) (see e.g. [1]) - or equivalently the normalised nuclear gluon distribution ratio 
R a (x,Q 2 ) = G (x,Q 2 )/A G P (x,Q 2 ) - is determined with a poor accuracy, say ~ 20%, 
and only for rather large values of x, x > 10~ 2 . As a consequence, the precision level for 
perturbative QCD cross sections in nuclear collisions (p-A or A-A) presently taken at 
RHIC or soon to be taken at the LHC is far from that achieved in p-p and p-p collisions 
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at Tevatron and LHC Let us also mention that probing accurately G A (x,Q 2 ) is also 
interesting in itself in order to shed light on non-linear evolution expected in QCD at 
small x; a field which has recently received a lot of attention (see e.g. [2]). For these 
reasons, any way leading to a better control of the gluon nuclear density needs to be 
considered and properly quantified. Based on a recent study [3], we explored in this 
talk how the production of prompt photons in hadron-nucleus reactions at RHIC and 
LHC may serve as a sensitive probe of R G (x, Q 2 ). 

At leading-order (LO) in a s , large-p x prompt photon production proceeds either 
via the annihilation process, qq — > gj, or by Compton scattering, q(q)g — > q(q)j- Note 
however that, below x ± = 2p ± /y / i < 0.1, the former channel is less than 10% of the 
latter. For clarity purposes, let us consider it to be negligible in the following. In the 
collinear factorisation framework, the inclusive production cross section in p-A collisions 
can simply be expressed as a product of parton distribution functional (PDF), in the 
proton and in the nucleus, times the partonic cross section, a: [3] 



dy d 2 p A 



dv 

■: ± eV/2 
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The sum over the quark flavours is implicit in the definition of F(x) = F2(x)/x. 
The integral over v, stemming from the rapidity of the unobserved gluon recoiling 
jet, y jet =asinh[t>~ 1 — (1 — therefore spoils the simple relationship between the 

observable on the one hand, and the nuclear parton densities - which we would like to 
constrain - on the other hand. However, as noticed in [3], the fairly slow variation of 
R G (x) and R p (x) with x as compared to that of the PDFs themselves allows Eq. (TTj) 
to be approximated by: 

dV(p A - 7 X ) _ , rA r vpP f x_^\ QP f xj_ 



dy d 2 p ± GV ± 'J V 2v J V 2 (l - v 

-,A / —7i\ f ^ s-yP ( % I ^ \ t-tP ( 30 , C 



= ^(^e-*) I„(y) + <(i ± e-») I„(y), 

where we used the fact that the integral is dominated by v ~ 1/2. Within this 
approximation, the inclusive prompt photon cross section is simply a linear combination 

A A 

of R and R p , both being evaluated at x ± e~ y . The above relationship simplifies 
even further in three kinematic limits. At mid-rapidity, I qg = I gq , while we have 
I qg (y) S> Igq(y) at forward rapidity, and I qg {y) Igq{y) a t backward rapidity. The 
last two inequalities arise from the dynamics of partons inside the proton: large-x 
is dominated by valence while gluons are dominant at small x. Therefore, photons 
produced at large positive (respectively, negative) rapidity necessarily originate from 

§ Parton distributions need to be evaluated at a factorisation scale O (p ± ) . This scale dependence does 
not play an important role in the present discussion and is therefore not made explicit in the following. 
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the scattering of a valence quark inside the proton (nucleus) and from a gluon inside 
the nucleus (proton). Defining the nuclear production ratio as 



it immediately follows that: 

• at mid-rapidity: R pA (x ± ,y — 0) ~ 0.5 R^(x ± ) + 0.5 R^(x ± ); 

• at forward rapidity: R pA (x ± e~ y ) ~ R^(x ± ); 

• at backward rapidity: R pA (x ± e~ y ) ~ R A p (x ± ). 

We identify at least two reasons which could somehow complicate the above 
relationships. First of all, photons may also come from the collinear fragmentation 
of hard quarks and gluons produced in the partonic process. In that case, the photon 
transverse momentum can no longer be related to the momentum fractions of quarks and 
gluons inside the incoming hadrons, and hence to the proton/nucleus parton densities. 
However, a way to get rid of most of these fragmentation photons is to apply isolation 
criteria. In that sense, the production of isolated photons is a cleaner probe of R Q and 
R p rather than the inclusive yield. Furthermore, the above approximations have been 
derived at leading order, which may well be spoiled when higher order corrections are 
taken into account. 



2. Phenomenology at RHIC 

In order to quantify the accuracy of the above approximations, the nuclear production 
ratio of isolated photons at forward rapidity (taking y — 3, a region which should be 
accessible both by the PHENIX [4] and the STAR [5] experiments at RHIC) is computed 
at NLO (using the calculations of Ref. [6]) in d-Au over p-p collisions at y/s NN = 200 GeV 
(Fig. ED. 

The expected suppression using the nDSg NLO [7] nuclear parton densities (solid 
line) is R pA ~ 0.65, that is much below R ~ 0.8 (dashed line). The main reason for such 
a mismatch is, however, coming from an isospin effect, as we compare systems which 
projectiles differ: a proton and a deuteron. Recall indeed that, being electromagnetic 
probes, photons couple more to u quarks than to d quarks (e^ = 4e^). Therefore, since 
the deuteron is poorer in u than the proton is, prompt photon production is suppressed 
in d-p with respect to p-p scattering. The proton and neutron seas being roughly 
symmetric, the suppression should be visible as long as valence quarks are involved, 
that is at large x\ > 0.1 which is the case at RHIC at forward rapidity. For the sake 
of a rough estimate, let us assume that only scattering of valence quarks with gluons is 
involved at large x and that u p ~ 2d p , then the isospin-corrected relationship we could 
expect between R dA and R Q is: R dA {x) ~ (5m + 5d)/(8u + 2d) x R (x) ~ 5/6 x R G (x) 
(in practice the isospin correction factor is given by R dp = a(dp)/a(pp), which is under 
better theoretical control). As can be seen in Fig. [1] (dotted line) the isospin-corrected 
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is in pretty good agreement, say within 10%, with the predicted R dA (see also the 
ratio of these two curves, minus one, r^ HIC , bottom line). 
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Figure 1. R dAu of isolated photon production Figure 2. Counting rates for isolated photons 
(y=3) at ^/s~^ = 200 GeV (see text for details) . in d-Au collisions at — 200 GeV at 

RHIC I and RHIC II. 

From this analysis [3], it appears that photon production is a sensitive probe of 
the nuclear modifications of parton densities. This observable is, nevertheless, only 
interesting if the nuclear production ratio can be measured with a high statistical 
accuracy, that is if reasonable rates can be achieved. Using an integrated luminosity of 
C = 0.02 pb" 1 and £ = 0.45 pb" 1 at RHIC I and RHIC II [8], respectively, the expected 
NLO rates are plotted in Fig. [2] as a function of p ± in the y = 2.5-3.5 rapidity range. 
Fortunately, the photon yield is pretty large up to rather high p ± and already at RHIC I, 
which is very encouraging. Provided the experimental systematic uncertainties remain 
low, there is a real hope that the suppression of isolated photon production soon to be 
measured at RHIC will allow for tight constraints to be set on G A /G P at rather small 
values of x. 
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